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Abstract 

In order to fully address the oscillation processes at a neutrino factory, a de- 
tector should be capable of identifying and measuring all three charged lepton 
flavors produced in charged current interactions and of measuring their charges 
to discriminate the incoming neutrino helicity. This is an experimentally chal- 
lenging task, given the required detector mass for long-baseline experiments. We 
address the benefit of a high-granularity, excellent-calorimetry non-magnetized 
target-detector, which provides a background- free identification of electron neu- 
trino charged current and a kinematical selection of tau neutrino charged current 
interactions. We assume that charge discrimination is only available for muons 
reaching an external magnetized-Fe spectrometer. This allows the clean classifi- 
cation of events into electron, right-sign muon, wrong-sign muon and no-lepton 
categories. In addition, high granularity permits a clean detection of quasi-elastic 
events, which by detecting the final state proton, provide a selection of the neu- 
trino electron helicity without the need of an electron charge measurement. From 
quantitative analyses of neutrino oscillation scenarios, we conclude that in many 
cases the discovery sensitivities and the measurements of the oscillation param- 
eters are dominated by the ability to measure the muon charge. However, we 
identify cases where identification of electron and tau samples contributes signif- 
icantly. 
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1 Introduction 



The firmly established disappearance of muon neutrinos of cosmic ray origin |T] strongly 
points toward the existence of neutrino oscillations || . The first generation long base- 
line (LBL) experiments - - K2K [§, MINOS @, OPERA § and ICANOE [§ - 
will give a conclusive and unambiguous signature of the oscillation mechanism and 
will provide the first precise measurements of the parameters governing the oscillation 
mechanism. MiniBOONE and the LBL programs will test the LSND signal ||. 



A neutrino "factory" [[], |TD[ is based on the decay of muons circulating in a storage 
ring. Neutrino factories raised the interest of the physics community, since they appear 
natural follow-ups to the current experimental LBL program and could open the way 
to future muon colliders. 



As many studies have shown [11], [12], [L3| 14]], the physics potential of such facili 



ties are indeed very vast. An entry-level neutrino factory could test the LSND signal 



in a background free environment ||12|| . More importantly, a neutrino factory source 
would be of sufficiently high intensity to perform very long baseline (transcontinen- 
tal) experiments. It could also bring the neutrino sector into the realm of precision 
measurements. 

The neutrino oscillation phenomenology may be complicated and involve a combi- 
nation of transitions to z/ e , and v T . It is quite evident that future neutrino factories 
will provide ideal conditions for the neutrino oscillation physics |fLT| , |13| , [14] , I5[ |I6 |. The 



neutrino flavor phenomenology could be completely explored: a precise measurement 
of the mass difference and mixing matrix elements is achievable, a test of the unitarity 
of the mixing matrix can be performed, a direct detection of Earth matter effects is 
feasible and CP violation effects could be studied on the leptonic sector |T5| ]. 



The combination of data from atmospheric neutrino and first generation LBL exper- 
iments will provide some preliminary information on the possible sub-leading electron 
mixingQ. A neutrino factory can largely improve the sensitivity on this mixing angle. 

Neutrino sources from muon decays provide clear advantages over neutrino beams 
from pion decays. The exact neutrino helicity composition is a fundamental tool to 
study neutrino oscillations. It can be easily selected, since /i + — > e + u e u^ and pT — > 
e'VeiSp can be separately obtained. 

At a neutrino factory, one could independently study the following flavor transitions: 
H~ -» e~ v e 

— > v e — > e~ appearance (1) 

— > disappearance, same sign muons (2) 

— ► v T — > appearance, high energy nu's (3) 

— > v e disappearance (4) 

— > — ► /i + appearance, wrong sign muons (5) 
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— > v T — ► r + appearance, high energy mi's (6) 
plus 6 other charge conjugate processes initiated from /i + decays. 

The other main advantages over traditional pion beams are (1) the beam is free 
of systematics and the composition is well known, therefore ideal for disappearance 
studies; (2) the two neutrinos in the beam have opposite helicities, therefore one can 
envisage oscillation appearance searches without intrinsic beam backgrounds (3) the 
muon energy is monochromatic and in principle adaptable (4) muon storage rings 
allow for multiple baselines, and hence a complete exploration of the L/E domains 
of oscillations (5) because very high intensity will be needed for muon colliders, very 
intense muon sources will produce very intense neutrino sources, at least a factor 100 
more intense than existing high energy facilities. 

While physics motivations are well understood, it is not yet clear which design of 
detector would best allow to take full advantage of the neutrino factory beams. 

We think that, in order to fully explore the neutrino oscillation processes, the 
detector should be capable of: 

1. measuring and identifying all three lepton flavors: electron, muon and tau; 

2. measuring the sign of the lepton charge; 

3. separating between charged and neutral current interactions. 

Experimentally, it is a very challenging task to build detectors with (1) mass scales 
of the order of tens of ktons required for long-baseline experiments, (2) with sufficient 
granularity to cleanly identify electron and tau leptons and (3) which measure the 
charge of these leptons. 

Various solutions have been explored recently. One based on nuclear emulsions and 
magnetized iron has been discussed in [[17]]. The main challenge there is to reach the 



required mass. Large magnetized calorimeters have been discussed in [|T^]. Such high- 
density detectors, while "easily" conceived as massive objects, have intrinsically very 
coarse granularity and only allow the clean measurement of muons. They certainly do 
not have sufficient power to adequately identify and measure electron or tau charged 
current states. 

In this paper, we are motivated by the recent progress made in the direction of 
the design of the multikton ICANOE detector 0: accordingly, we consider a 10 kton 
(fiducial) high granularity low density liquid argon imaging target, complemented with 
a high-acceptance external muon spectrometer. 

Thanks to its extremely high granularity target and its excellent calorimetric prop- 
erties, this design provides the clean identification and measurement of all three neu- 
trino flavors: electron, muons and taus. However, only the sign of the muons reaching 
the muon spectrometer can be determined^. 

1 The possibility of the measurement of the electron charge will be addressed in a future work. 
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The aim of this paper is to understand the potentials of a non-magnetized high- 
granularity target detector which, compared to traditional high density iron calorime- 
ters, brings the measurement of electrons and taus. 

We study the physics potentials of such a detector configuration for three-family 
neutrino mixing and for three different baselines (732, 2900 and 7400 km). 

In Section we summarize three-family neutrino mixing framework, including the 
treatment of propagation through matter. 

In Section |3], we explain in details how the event distributions and rates are obtained 
from a detailed simulation of neutrino interactions and detector effects. 

In Section ^, we construct for given oscillation scenarios, event variable distributions 
for various event classes. All events can be subdivided into the electron, the same sign 
muons, the wrong sign muons and the no lepton samples. 

In addition, in Section |5] we discuss the possibility to further discriminate final states 
between v e and origins from v T by means of kinematical analysis of the events. 

We also address in Section || the possibility to tag quasi-elastic events which provide 
indirect neutrino helicity discrimination however at a large statistical price. 

Section [7] is devoted to describing our fits of the oscillation parameters. The fits 
are expected to give back the input reference oscillation parameters and are used to 
estimate the precision with which we can estimate these parameters. 

Section |8] presents the results for the important case in which the oscillation effects 
can be approximated by one mass scale. 

Since the information about the oscillation parameter is redundantly available in 
the visible energy distributions of the various event classes, we address in Section || the 
question of the consistency between the different observed oscillations processes. The 
ability to treat the appearance of electron or tau neutrinos gives good over-constraints 
on the mixing matrix. 



Finally, in Section [10], we analyze the three-family scenario including possible CP- 
violation. 



2 Three-family neutrino oscillation framework 
2.1 Mixing matrix parameterization 

We consider neutrino oscillations in a three-family scenario: the flavor eigenstates 
v a {a = e, fi, r) are related to the mass eigenstates u[{i = 1, 2, 3) by the mixing matrix 
U 

v a = XJ ai v[ (7) 
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and we parameterize it as: 

/ C12C13 S12C13 s 13 e~ lS \ 

U = -S12C23 - Ci 2 Si 3 S23e i ' 5 C12C23 - Si2Si3S2 3 e l<5 C13S23 (8) 

I iS iS I 

\ S12S23 _ c 12Sl3 c 23 e — C12S23 — S 12 Si3C23e C13C23 / 

with S{j = sin 9^ and qj = cos 9^. We confine without loss of generality the mixing 
angle 6*13 to values in the interval [0, 7r/4], and 6*12, 6*23, 5 to the interval [0, 7r/2]. We 
present results in terms of sin 2 29 13, sin 2 823 and sin 2 # 12 , all running in the interval 
[0,1]. The reason for this choice can be for example seen in Appendix A, where we 
recall oscillation probabilities in the one mass scale approximation. The oscillation 
probabilities for v e — > and v e — > v T depend on sin 2 29 \$ and on the sin and cosine 
of #23- The angle #23 must span the interval [0,7r/2], however, the #13 can vary within 
interval [0, 7r/4]. Note that we consider small 6*13 angles, so the cos 4 6*13 dependence in 
V/i v r is very mild. 

For 5 = (i.e. U is real), the general expression for the three-family neutrino 
oscillation probability is: 

P{y a -> v P ] E, L) = P(v a -> vp\ E, L) = 5 af3 - 4 ^ J af3jk sin 2 (A jk ) (9) 

3>k 

where in natural units = Am 2 fc L/4i? = (m 2 — m1)L/4E, E is the neutrino energy, 
L is the neutrino path- length, and the Jarlskog term is J a /3jk = UpjUpkU a jU a k- 

We naturally assign the mass difference squared Am 2 2 to explain the solar neutrino 
deficit and the mass difference squared Am 2 2 to describe the atmospheric neutrino 
anomaly. We will take as reference value Am| 2 = 3.5 x 10 -3 eV 2 . To cover possible 
ranges of this value, we will also consider two other values Am| 2 = 5 x 10 -3 eV 2 and 
Am| 2 = 7 x 10~ 3 eV 2 . We will always assume maximal (2-3)-mixing sin 2 2# 2 3 = 1- 
In case of the solar neutrino deficit solution, the values for Am 2 2 and mixing sin 2 29 12 
are not uniquely defined by experiments. We will limit ourself to the LMA-MSW 
solution with parameters Am 2 2 = 1 x 10~ 4 eV 2 and hypothesize a maximal (1-2)- 
mixing sin 2 9 12 = 0.5. 

In this paper, we do not consider the LSND result which would force us to include 
more states with new parameters beyond three-family mixing. 

This choice implies that we will always work in a situation where | Am 21 1 < | Am 2 2 1 w 
jAmfiJ. In first approximation, the oscillation phenomena governed by the two mass 
differences decouple and the effects produced by Am 2 2 are small at high energy for the 
considered baselines. In the first part of this paper, we will neglect Am 2 2 effects and 
work in the so-called "one mass scale approximation" [PI . In a second phase, we will 
be concerned with CP- violation effects and will have to include Am 2 2 . 

For simplicity, we will take m\ < m 2 < 772,3 which implies Am 2 2 > 0. We recall 
that neutrino oscillations through matter can be used to distinguish Am| 2 > from 
AmL < 0. 
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In vacuum, we will express the oscillation probability as a function of the seven fol- 
lowing parameters: (a) the three mixing angles #12, #13, #23; (b) the two mass differences 
squared Am 2 2 , Am^; (c) the baseline L; (d) the neutrino energy E. 



2.2 Matter effects 

Since we will consider very long distances between neutrino production and detection, 
this will only be possible in practice for neutrinos traveling inside the Earth. In this 
case, the neutrino oscillation probabilities will be modified by an additional diagram 
due to the interaction of electron neutrino with the electrons in the matter PH. One 
can maintain the neutrino oscillation formalism derived in vacuum but define effective 
masses and mixing angles valid in matter. For example, the effective masses will result 
from the diagonalization of the Hamiltonian: 



U 

where 
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D = 2V2G F n e E = 7.56 x lO" 5 ^— ?—)(——) (11) 

gem 6 Gel/ 

Here, n e is the electron density and p the matter density. For anti-neutrinos, we must 
replace D by —D. For p(g cm~~ 3 )E(GeV) ~ 40, the effective mass parameter D is of 
the order of the mass splitting that is derived from the atmospheric neutrino anomaly. 
We then expect matter effects to be important. 

Within the two-family mixing scheme, the modification of the flavor transition in 
matter is taken into account by the mixing angle in matter 9 m , which is: 

9 „ / x sin 2 29 , . 

sin 2 29 m (D) = . 2 5 — 2 12 

sm z 29 + {-^ - cos 29) 1 

For neutrinos, a resonance condition will be met when D ~ Am 2 cos 29 and the oscil- 
lation amplitude will reach a maximum. The resonant neutrino energy E res is 

^ 1.32 x 10 4 cos 29Am 2 (eV 2 ) 

h n — 3^ ( 16 ) 

Rather than two-family mixing, we have adopted throughout this study three-family 
framework. In this context, we use the analytic expressions for the matter mixing angles 
and mass eigenvalues calculated in |23], |2l| (see Appendix A). 

The mass eigenvalues in matter Mi, M2 and M3 are: 

Ml = m \ + -- -VA 2 - 3BS - —VA 2 - WVT^S 2 (14) 
3 3 3 
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Figure 1: Evolution of mixing angles in matter as a function of D(= 
2y/2G F N e E u )/Am 2 2 . The reference vacuum parameters are given in the figure. 



Ml = ml + j - ^VA 2 - 3BS + ^j-VA 2 - wVl^S 2 (15) 
A 



Ml = m? + - + -VA 2 - 3BS (16) 
6 1 3 3 

where A, I? and S are given in the Appendix. For the mixing angles in matter the 
analytical expressions read: 

• 2 Q m -(M^aM 2 + P)AM 2 1 

12 ' AM| 2 (Mf - aM 2 + (3) - AMfj (M 2 4 - aMf + /3) 1 J 

Sm ^ " aMam| (18) 

• 2fl m G 2 g| 3 + ^4, + 2GFc 23g23 c, 
sm c/ 23 - Q2 + p2 

where a, (3, G and F are found in the Appendix. 
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To illustrate matter effects in three-neutrino mixing framework, we show in Figure |I| 
the values of the mixing angles in matter, plotted as a function of D/Am^, or equiv- 
alently of p x E. The parameter values in vacuum correspond to our reference values 
for atmospheric and LMA-MSW solar experiments. The resonant behavior of 9™^ is 
clearly visible. It gives maximum oscillation at a neutrino energy of about 12 GeV for 
a density of 3.7 g cm~ 3 . There is a similar resonant behavior for 8^ but it occurs at 
low energy since it is driven by Am^. For D > Am| 2 , the angles sin 2 29^ and sin 2 6^ 
tend to rise slightly, because the non- vanishing Am^ splitting removes the degeneracy 
between muon and tau flavors. 

Our results have been computed assuming a constant density along the whole neu- 
trino path, and equal to the mean density, obtained integrating over the earth profile 



2jJ . This approximation yields, as shown in e.g. Ref. p5fl , similar results to those 
obtained by numerical integration using the actual Earth's density profile. 

The oscillation probability through matter will be a function of eight parameters: 
(a) the vacuum three mixing angles #12, #13, #23; (b) the vacuum two mass differences 
squared Am 2 2 , Am 2 3 ; (c) the average earth density p; (d) the baseline L; (e) the 
neutrino energy E. 



3 Choice of baseline and event rates 

The exact parameters of a neutrino factory are not yet completely fixed and realistic 
scenarios are in the process to be defined flOU- However, based on [|l(J, we can assume 
that the muons in the storage ring have an energy E^ = 30 GeV and that after one 
year of operation, the factory should deliver about 10 20 "useful" muons decays of both 
polarities in the straight section pointing towards the far detector location. We base 
our ultimate reach on 10 21 "useful" muons decays. Even an integrated intensity of 10 22 
might be eventually reachable. 

We compute the fluxes assuming unpolarized muons and disregarding muon beams 
divergences within the storage ring. We integrate the expected event rates using a 



neutrino- nucleon Monte-Carlo generator p8[. The total charged current (CC) cross 



section is technically subdivided into three parts: the exclusive quasi-elastic scattering 
channel oqe and the inelastic cross section Oi ne iasic which includes all other processes 
except charm production which is included separately. 

Table [I] summarizes the expected rates for the 10 kton fiducial mass and 10 20 muon 
decays (expected 1 year of operation). N to t is the total number of events and N qe is 
the number of quasi-elastic events. 

Even though our study is site non-specific, the chosen baselines could correspond to 
the distances between the Laboratori Nazionale del Gran Sasso (LNGS) and neutrino 
factories at (1) CERN (L = 732 km, < p Ear th >= 2.8 g/cm 3 ), (2) Canary Islands (L = 
2900 km, < pEarth >= 3.2 g/cm 3 ) and (3) Fermilab (L = 7400 km, < p Earth >= 3.7 
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Event rates for various baselines 
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i/ e cc 


197000 


7880 


12900 
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z/ e NC 


57900 




3670 




580 





Table 1: Expected events rates for a 10 kton (fiducial) detector in case no oscillations 
occur for 10 20 muon decays. N tot is the total number of events and N qe is the number 
of quasi-elastic events. 



g/cm 3 ). 

4 The four main classes of events 

Muon identification, charge and momentum measurement provide discrimination be- 
tween Up and charged current (CC) events. Good v e CC versus v NC discrimination 
relies on the fine granularity of the target. Finally, the identification of v T CC events 
requires a precise measurement of all final state particles. 

It is natural to classify the events in four classes |TTJ. We illustrate them for the 
case of fi~ stored in the ring. 

1. Right sign muons (rs/x): the leading muon has the same charge as those 
circulating inside the ring. Their origin is from 

(a) non-oscillated CC 

(b) — > v T CC, r~ — > fi~ decays 

(c) hadron decays in neutral currents. 

2. Wrong sign muons (wsn): the leading muon has opposite charge to those 
circulating inside the ring. Opposite-sign leading muons can only be produced 
by neutrino oscillations, since there is no component in the beam that could 
account for them. 

(a) u e — > oscillations 
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(b) v e — > v T oscillations, r + — - > /i + decays 

(c) hadron decays in neutral currents. 

3. Electrons (e): events with a prompt electron and no primary muon identified. 
Events with leading electron or positron are produced by the charged-current 
interactions of the following neutrinos: 

(a) non-oscillated z7 e neutrinos 

(b) — > v e oscillations 

(c) v e — > v T or Up — > v T oscillations with r — > e decays 

4. No Lepton (0£): events corresponding to NC interactions or v T CC events 
followed by a hadronic decay of the tau lepton. Events with no leading electrons 
or muons will be used to study the — > v T oscillations. These events can be 
produced in 

(a) neutral current processes 

(b) v e — > v T or — > v T oscillations with r — > hadrons decays 

The last two classes can only be cleanly studied in a fine granularity detector. 

The most effective way to fit the oscillation parameters is to study the visible energy 
distribution of the four classes of events defined above, since assuming the unoscillated 
spectra are known, they contain direct information on the oscillation probabilities. 

Of course, for electron or muon charged current events, the visible energy recon- 
structs the incoming neutrino energy. In the case of neutral currents or the charged 
current of tau neutrinos, the visible energy is less than the visible energy because of 
undetected neutrinos in the final state. The information is in this case degraded but 
can still be used. 



Our analyses are performed on samples of fully generated Monte-Carlo |28|] events, 
which include proper kinematics of the events, full hadronization of the recoiling jet 
and proper exclusive polarized tau decays when relevant^. Nuclear effects, which are 



taken into account by the FLUKA model |[29|| , are included as they are important for 
a proper estimation of the tau kinematical identification. 

The detector response is included in our analyses using a fast simulation which 
parameterizes the momentum and angular resolution of the emerging particles, using 
essentially the following values: electromagnetic shower 3%/v^E©l%, hadronic shower 
~ 20%/\/E © 5%, and magnetic muon momentum measurement 20%. 

Hadron decay background can be quite large in a low density target and could 
be quite dangerous. Fortunately, it can be easily suppressed by a cut on the muon 



2 Our simulation has been bench-marked on the comparison of the kinematic features of the lepton 
and hadronic jet of real neutrino data accumulated in the NOMAD experiment (see e.g. ||). 
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candidate momentum, > 2 GeV, which reduces it to a tolerable level. Figure § 
illustrates the relative background expected for and NC and CC processes as a 
function of the muon momentum. After the cut, the expected contamination for 
CC events is at the level of 10~ 5 . Real charged current events maintain an efficiency 
above 95%. 

The visible energy is computed as the modulus of the vector sum of the momenta 
of each visible particle in the event. Figures [| |], |5] and |S| show the reconstructed 
visible energy at the baseline L = 7400km normalized to 10 20 //'s for each event class 
for a specific oscillation scenario with Am 2 2 = 3.5 x 10~ 3 eV 2 , sin 2 9 2 3 = 0.5 and 
sin 2 26 , 13 = 0.05. The different contributions including backgrounds for each event class 
have been evidenced in the plots. For example, in Figure |], the different processes that 
contribute to the right-sign muon class are unoscillated muons, taus and background 
events. 



5 Further classification based on kinematical anal- 
ysis 

An efficient identification of v T induced charged current events requires a precise mea- 
surement of all final state particles. Excellent calorimetry allows to take full advantage 
of the special kinematic features of v T events. We independently search for the leptonic 
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Figure 3: Visible energy spectrum for elec- 
tron events: v e CC (dashed line), v T and 
v T (dotted line) and v e CC (dot-dashed). 
The solid histogram shows the sum of all 
contributions. 
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Figure 5: same as Figure |] for wrong 
sign muon sample: CC (dashed line), 
v r and z/ T (dotted line) and meson decay 
background (dot-dashed). The solid his- 
togram shows the sum of all contributions. 



Figure 4: same as Figure |^ for right-sign 
muon sample: CC (dashed line) , v T and 
v T (dotted line) and meson decay back- 
ground (dot-dashed). The solid histogram 
shows the sum of all contributions. 
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Figure 6: same as Figure ^ for the no- 
lepton sample: v NC (dashed line), v T 
and v T (dotted line). The solid histogram 
shows the sum of all contributions. 
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v T appearance search 


Cuts 


r -> / 


CC background 


Cuts 




NC background 


Initial 


100% 


100% 


Initial 


100% 


100% 


Loose cuts 


P T < 0.5GeV 


50% 


14% 


pmtss < 1Ge y 


72% 


40% 


P™ iss > 0.6GeV 


40% 


0.5% 


Q T > 0.5GeV 


30% 


2% 


Tight cuts 


P T < 0.5GeV 


50% 


14% 


pmtss < 1Ge y 


72% 


40% 


pmiss > 1Ge y 


20% 


0.08% 


Q T > lGeV 


6% 


0.07% 



Table 2: z/ T appearance search for the leptonic and hadronic decay modes of the tau lep- 
ton. Overall signal efficiencies and fractional remaining backgrounds are quoted. The 
case labelled as "Tight cuts" correspond to the situation where one event background 
is expected, for 10 20 muon decays, at the farthest location (L = 7400 km). 



and hadronic tau decay modes. 

For the r — > Ivv decay mode, the main background comes from v% CC. To enhance 
the separation between r and background events, we demand the event missing P? to 
be larger than 0.6 GeV and the transverse momentum of the lepton candidate, P T , to 
be smaller than 0.5 GeV. This set of cuts is referred to as "loose cuts" and it will be 
used to perform a check on appearance/disappearance consistency (see subsection [5]). 
In Table 0, we show that the overall r efficiency for "loose cuts" is 40% for a CC 
background level of ~ 5 x 10~ 3 . Figure [7| shows the energy spectra for the four event 
classes after application of these cuts, for u T CC and other types of events. No energy 
cut has been applied, but the fact of using the energy spectra in the fit also exploits 
the difference in energy spectra of v T CC events. 

A set of "tight cuts" is also applied, aiming at having one expected background 
event at the farthest location in case 10 20 "useful" muons decays are delivered. As we 
can see from Table |2], the overall tau efficiency in this case amounts up to 20%. 

For hadronic decays the most important source of background correspond to NC 
events. If we demand a P™ tss smaller than 1 GeV and a transverse momentum of the 
hadron candidate with respect to the total event momentum, Qt, larger than 0.5 GeV 
("loose cuts") only 2% of the initial background survives for an overall tau efficiency 
of 30%. If we require only one NC background event survivor at L — 7400 km ("tight 
cuts"), the signal efficiency drops to 6% due to the stringent Qt requirement imposed. 
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Figure 7: Visible energy spectrum for the four event classes after application of loose 
kinematic cuts (see text). The presence of r events in the third class (wrong-sign 
muons) is an indication of the process v e —>■ v T . 
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6 Quasi-elastic final states 



The quasi-elastic process, while rare, is a clean process that allows to separate neutrino 
from anti-neutrino events, in principle for all neutrino flavors, since Ug + n — > + p 
and i>i + p — > £ + + n. The recoil proton is easily identifiable within the high-granularity 
target. 

This channel is particularly interesting to study oscillation in the electron channel. 
Starting from negative muons circulating in the storage ring, we look for exclusive 
electron-proton final states. These provide "background-free" oscillation signals: 

-> v e + n -> e" + p (20) 
^v T -\-n^T~+p—>- e~vv + p (21) 

since v e + p — > e + + n. 

The selection of quasi-elastic events is the only way to identify the helicity of neu- 
trino electrons in absence of measurement of the electron charge. 



7 Oscillation parameters fitting 

Given the adopted parameterization of the mixing matrix, we have a priori a total of 
7 free parameters, which can be represented by the vector: 

P = (Am 2 1; Am 2 32 , sin 2 12 , sin 2 20 13 , sin 2 # 23 , 5, p) (22) 

The values of the parameters governing the oscillations are extracted from a global 
fit of the visible energy distributions obtained for each event class. The fit is performed 
with the MINUIT [22| package, and is expected to get back the same values of the 



parameters, starting from the reference distributions. 

At each iteration, a different set of parameters is probed, and with the same pro- 
cedure used to get the reference histograms. 

For a given polarity A of the muons in the storage ring, we compute x 2 's of the 
difference between the binned oscillated spectra, which will be function of the parame- 
ters, and the reference histograms. We define a x 2 f° r each of the four classes of events, 
i.e. the electrons (e), the right-sign muon (rs^), the wrong sign muons (tus/i) and the 
no lepton class (0£): 

+ xl,oi (23) 



where 



N9(P) - N[(P ref) 



2 



XUP) = E[ tV ' ; wre " l • (24) 

0*« r 
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The sum runs over 25 equally spaced energy bins (E^ i — 1, 25). Here, N?(P re f) is the 
expectation for the reference values and AT?(P) is the "data" obtained for a given set of 
the oscillation parameters P. o i c contains both statistical and systematic contributions 
to the total error. In order to improve the statistical treatment of bins with low statis- 
tics, a bin that possesses less than 40 events is assumed to be Poisson distributed and 
therefore its contribution is computed as 2(Nf(P)-Nf(P ref ))+2N i \n{Nf{P) / N?{P ref )) 
(see Ref. §(J). 



The systematic error takes into account the uncertainties in the knowledge of the 
beam, neutrino cross sections and selection efficiencies and we assume it amounts up 
to 2% uncorrelated from bin to bin. 

We assume that a neutrino factory will operate with alternate runs of opposite 
muon polarities, therefore eight energy distributions can be fitted simultaneously: 

Xall — X+,all + X-,all (25) 

The values of the fitted parameters P are obtained minimizing the \ 2 {P)- 

It is in practice not always possible to fit all the free parameters, since for some 
parameter-space regions, the oscillation effects at the chosen baselines and energies 
can be negligible. In particular, this can be the case for Am^ and sin 2 9\2 which drive 
the solar oscillations. For values Am^ <C 1CT 4 eV 2 , we are insensitive to the "solar" 
sector. 

We therefore adopted successive fitting procedures with an increasing number of 
free parameters. 

At first, we can simplify the three- family oscillation picture if the oscillations pro- 
duced by Am^ can be neglected at the considered baselines and energies. In this 
case, only the mass difference squared Am^ and the two mixing angles #13 and #23 
are relevant. The oscillation probabilities are given in the Appendix. For example, for 
—>■ Vy, oscillations, it is: 

P{v e - V(t , E, L) = sin 2 (2^ 3 ) sin 2 (fl™) A 2 2 (26) 

where A§ 2 = sin 2 ((Mf - M 2 2 )L/4£). The fit has in this case 4 free parameters, which 
can be represented by the vector: 

P lms = (Am 2 2 , sin 2 20 13 , sin 2 6 23 , p) (27) 

The S phase has disappeared, since within this approximation it becomes unphysical. 
The results are presented in Section §, assuming a Am 2 2 parameter varying in the range 
of values favoured by current atmospheric data (Am 2 2 = 3.5, 5, 7x 10~ 3 eV 2 ), a maximal 



(2-3)-mixing sin 623 = 0.5 and a #13 value compatible with CHOOZ results [26 and 
recent fits to data @ (sin 2 2# 13 = 0.05). 

We return to the general three-family scenario in Section |1(] where we consider 
sensitivity to the CP violation. 
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Two family mixing 

5(sm 2 6 23 ) for # 23 = 45°, 13 = 




Only right-sign muons 

X,±,rsfi 


All classes 

X.±,all 


Am| 2 (eV 2 ) 


L=7400 km 


L=2900 km 


L=7400 km 


L=2900 km 


7 x 10~ 3 
5 x lO- 3 
3.5 x 1(T 3 


22% 
12% 
10% 


8% 

12% 

18% 


20% 
11% 
10% 


8% 

10% 

16% 



Table 3: Precision in the measurement of the mixing angle assuming two family mixing 
for three possible mass differences and two very large baselines. In all the cases the 
precision obtained in the measurement of Am^ is 1%. 



In order to compute the precision of the determination of the parameters, we con- 
sider two methods: (1) a one-dimensional "scan" of a given parameter; the other 
variables are left free and mininized at each step; the resp. 1,2,3 sigmas are given by 
resp. Xmin + 1) +4 an d +9. (2) a two-dimensional "scan" of a two-parameter plane; 
the other variables are left free and minimized at each point in the plane; the resp. 
68%, 90%, 99% C.L. are given by resp. xLn + 2 -3, +6.0 and +9.2. 



8 Results for one mass scale approximation 
8.1 Case of two- family mixing : #13 = 

The detection of the dip in the energy distribution of the right-sign muon sample 
dominates the precision on the measurement of the mixing angle #23 and the mass 
difference Am^. This is true provided that the beam energy and baseline are chosen 
in such a way that the disappearance maximum is visible in the oscillated spectrum. 
Table |3| summarizes the expected accuracies. With 10 20 muon decays of each polarity, 
precisions of 1 — 2% are expected in the determination of Am| 2 while for the precision 
on sin 2 #23 is around 10%, in agreement with results quoted in Jl4 . 



We compare in figure |] how the precision on the mass difference and mixing angle 
changes when experimental resolutions and backgrounds are disregarded. We see that 
our fits are barely affected and only for the cases where the dip is not seen, the instru- 
mental effects and backgrounds spoil at the level of a few per cent the accuracy on the 
oscillation parameters. 
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Figure 8: Two-dimensional 90%C.L. contours for Am| 2 and sin 2 623 for three differ- 
ent Awigj values. Inner contours: no background and perfect muon resolution; outer 
contours: backgrounds and muon resolution included. 



8.2 Case of three- family mixing : #13 7^ 

In order to obtain the best sensitivity on the mixing angle #13, the search for wrong 
sign muons is ideal, since it will be a direct signature for v e — ► oscillations. 

For this kind of study, since we are dealing with the smallest number of signal events, 
the sensitivity does strongly depend on the ability of rejecting background. Figure § 
shows for L = 7400 km, the sensitivity on #13 for two different muon normalizations 
(10 20 and 10 21 muon decays of each polarity). For each pair of values (Am| 2 , 6*13), the 
fit was performed leaving 623 free. 

In the obtaining previous plot, we did not on purpose apply strong background 
cuts, since we believe that very high rejection powers obtained on paper may not 
stand the proof of real experimental conditions, with non-gaussian behaviours, tails of 
distributions etc. 

To also give the maximum of sensitivity that can be obtained in principle, we also 
illustrated in figure ^] the effect that a background free environment would have in the 
expected sensitivity. At 90% C.L., we obtain sin 2 26 ) i 3 < 10 -3 — 10~ 4 depending on 
the number of muon decays. This represents two orders of magnitude improvement 
with respect to quoted sensitivities at CNGS. In case we consider 10 21 muon decays 
and backgrounds are reduced to a negligible level, the obtained sensitivity is consistent 
with the one quoted in 



For a background-free environment gives the best sensitivity, the amount of infor- 
mation added by other event classes (i.e. the electrons) is negligible. 
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Figure 9: Sensitivity on 6*i 3 
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Three-family mixing 








sin 2 9 23 


= 0.5, sin 2 2# 13 = 0.05 








All classes 


Only muons 




X,all 








L=2900 km 


L=7400 km 


L=2900 km 


L= 


7400 km 




Am 2 2 = 3.5 x 10~ E 


eV 2 , 






5(sin 2 # 23 ) 


1.4% 


0.9% 


1.4% 




0.9% 


14% 


8% 


16% 




9% 


5(sin 2 2# 13 ) 


15% 


10% 


17% 




15% 




Am 2 32 = 5 x 1(T 3 


eV 2 






5 (Am 2 ) 
5(sin 2 # 23 ) 


0.4% 
11% 


0.8% 

8% 


0.4% 
10% 




0.8% 
12% 


5(sin 2 2#i 3 ) 


11% 


9% 


14% 




16% 




Am 2 2 = 7 x 10~ 3 


eV 2 






S(Aml 2 ) 
S(sm 2 9 23 ) 


0.4% 


0.6% 


0.4% 




0.6% 


7% 


8% 


8% 




18% 


5(sin 2 2# 13 ) 


8% 


6% 


9% 




20% 



Table 4: Precision on the measurement of the oscillation parameters. 



8.2.1 Determination of # 13 

The measurement of #i 3 can profit from long baselines, since matter effects will enhance 
the oscillation signal. In presence of backgrounds, it is more favorable to enhance 
neutrinos signal even at the cost of the suppression of the anti-neutrino oscillations. 

Table |4] summarizes the expected precision on the measurement of the oscillation 
parameters. 

In this case, since for the chosen value of sm 2 2#i 3 = 0.05 the number of signal events 
is quite large, there is not any more a strong need for a background-free environment. 
Therefore, the inclusion of other event classes, like the electrons, can help to constrain 
the oscillation parameters. 

8.3 Sensitivity to #13 with quasi-elastic events 

An exclusive way of detecting the effects of a non vanishing #i 3 is through the ap- 
pearance of wrong sign electrons. Although, with the assumed detector configuration, 
there is no ability to directly measure the charge of the leading electrons, there is the 
possibility of disentangling final state electrons from positrons through the use of quasi- 
elastic events. We look for electron-proton final states. For example, in the target of 
ICANOE, a proton can be resolved if its kinetic energy is larger than about 100 MeV 



19 



u e appearance search with quasi-elastic 
Electron Class: Events for 10 21 p~ decays 


Baseline 


v e CC 


Up -> u e CC 


— ► v T CC, r — > e 


Total 


Elastic 


Total 


Elastic 


Total 


Elastic 


L = 732 km 
L = 2900 km 
L = 7400 km 


860000 
54300 
8300 


43000 
2700 
410 


2090 
1720 
960 


84 
70 
40 


3990 
3300 
1450 


110 

90 

40 



Table 5: Expected number of electron type events for 10 21 p~ decays. The three 
contributions to the total number of electron events are shown separately. Rates have 
been computed assuming oscillations with Am 2 2 = 3.5 x 10~ 3 eV 2 , sin 2 #23 = 0.5 and 
sin 2 20 n = 0.05. 



corresponding to a range of more than 2 cm. 

Table ||] shows the expected rates contributing to the electron class before any cut is 
applied. The background is twofold: (a) quasi-elastic v T CC events followed by r — ► e, 
(b) v e CC with an extra proton from nuclear origin. 

This last background can be estimated from data themselves studying the reaction 
VpP —>■ fi + n. Demanding a back to back electron-proton event topology with a proton 
kinetic energy in excess of 100 MeV, we estimate the expected background to be less 
than one event for an overall signal efficiency of 50%. 

The quasi-elastic channel provides, in a "background free" environment, 20 to 40 
gold-plated events depending on the selected baseline. It is very clean channel but is 
limited by statistics. 



8.4 Fit of the average Earth density parameter 

In matter, the amplitude of neutrino oscillation goes through a maximum for an energy 
given by equation [13|. Since 9\s is small, the MSW resonance peak is only a function of 



p and Awijj, This can be seen in Figure 10. Since the mass difference is constrained by 



the disappearance of right-sign muons, p is well-determined by the energy distribution 
of wrong sign muons. 

We extract the density from the fit, leaving it as a free parameter, as well as Am 2 2 , 
023 and 0i3. The precision on the determination of p depends on the baseline, as shown 
in table |], obtained considering 2 x 10 21 muon decays. For the longest baseline, where 
matter effect are large, a precision as good as 2% can be obtained. 



The influence of p in our fits is addressed in figure |TT| We see that for L = 7400 km 
and three different muon normalizations, the fact that p is either considered as a free 
parameter or fixed during the fit does not influence the accuracy in the determination 
of the mixing angles. 
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Distance (km) 


Density (g/cm 3 ) 


Relative error (%) 


732 


2.8 


18 


2900 


3.2 


10 


7400 


3.7 


2 



Table 6: Precision on the determination of p, from a global fit where also Am| 2 , 623 
and #13 are left as free parameters. This result has been obtained for 2 x 10 21 muon 
decays. 

9 Over-constraining the oscillation parameters 

The information about the oscillation parameter is redundantly available in the visible 
energy distributions of the various event classes. This allows us to address the question 
of the consistency between the different observed oscillations processes. 

Given the high statistical accuracy of the measurement, the consistency can be 
tested with good accuracy. In the three active neutrino mixing scheme, this implies 
that Yj y=e ^ T P{y x — > v y ) should be equal to one for x = e, /i, r and the same holds for 
anti-neutrinos. 

Other models can predict different values (i.e., oscillations to sterile neutrinos exist, 
the sum would be smaller than one). 

Let us concentrate on the oscillations into r neutrinos. In the case of negative 
muons in the ring, they can be originated from v^ — > v T or v e — > vt oscillations. The 
latter case is particularly interesting, since coupling a neutrino factory with a detector 
with t identification capabilities is probably the only way to identify and measure 
such a process. The v e — > vt can be revealed experimentally from the presence of r 
candidates in the wrong-sign muon sample, due to the opposite helicity of electron and 
muon neutrinos in the beam. 

To have a quantitative estimation of the consistency check of the various oscillation 
modes to r neutrino, we assign two global normalization factors, a and (3 to the oscil- 
lation probabilities Piv^ —>■ v T ) and P(v e — > v T ). If no new phenomena occur, these 
parameters should be exactly one. From a global fit to the visible energy distributions 
it is possible to extract the values of these parameters, and the precision obtainable on 
their measurement. 

To select r events from the background, still retaining a high efficiency on the signal, 
we apply a set of loose kinematic cuts (see table 0). 

Background levels of the order of one event can be reached applying tighter cuts, 
but in these cases the statistics is too small and the results obtained are slightly worse. 

Since the r lepton decays into muons, electrons and hadrons, we expect that a fit 
to all event classes would result in a remarkable improvement on the precision for the 
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Figure 10: Variation of the MSW reso- 
nance peak for wrong sign muons as a 
function of Earth's density. The plot is 
normalized to 10 21 p + decays. 
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Figure 11: 68%C.L. two-dimensional con- 
tours for sin 2 2^13 and sin 2 623- influence 
of p in the determination of the mixing 
angles for three different muon normaliza- 
tions and L = 7400 km. In the upper plot 
p is fixed during the fit, while in the lower 
one is taken free parameter. 



a parameter. Figure 12 shows how a determination improves as the different event 



classes are included in the fit. 

Table [7] shows the expected precisions in the determination of a and (3. We observe 
that for 10 20 decays, a is better determined (accuracy around 5%) for L = 2900km, 
however for 10 21 muons, the accuracy is about 1% regardless of the baseline and the 
mass difference and therefore oscillations into a sterile neutrino can be largely ruled 
out. 

The accuracy on /3, and therefore the first experimental evidence for v e — *> v T is 
much worst, given the smaller statistics available, since this oscillation probability is 
smaller than the corresponding — > v r by a factor sin 2 2# 13 , taken to be in this case 
0.05. We observe that somewhat better determination exists at L = 7400km since this 
oscillation mode is largely influenced by matter effects given its dependence on 13 . We 
conclude that a precision at the level of a few per cent on the observation of v e — > v T 
oscillations would require O(10 22 ) useful muon decays. 
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Figure 12: Precision on the determination of the a parameter when the fit is performed 
using: all event classes(solid line), right-sign muons + electrons(dashed line) and right- 
sign muons only(dotted line). 
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60% 


20% 
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7400 km 
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35% 


10% 


5% 
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25% 


7% 


2% 




3.5 


75% 


25% 


9% 


2900 km 


5 


25% 


15% 


5% 




7 


30% 


10% 


4% 



Table 7: Precision on the determination of the parameters a and (3 that quantify, 
respectively, the amount of — > v T and v e — » z/ T present in the data. In the three 
neutrino framework, the reference values are a, (3=1. 



10 Results for the general three-family scenario and 
CP-violation 

Let us consider now a more complex scenario. In this case, the value of the mass 
difference Atoj 2 is not any more negligible, and is actually assumed to be lO^eV" 2 , one 
of the highest values compatible with the large mixing angle MSW solution for solar 



neutrinos 21 



The oscillation does not depend any more on only three parameters, but all four 
independent angles of the mixing matrix and the two mass differences become impor- 
tant. 

In the most general case, the phase 5 can be different from zero, producing a complex 
mixing matrix, and thus generating CP violation. 

We recall that for neutrinos propagating in vacuum, the oscillation probability after 
a distance L for neutrinos can be expressed as: 

P{v a ^v p -E,V) = P CP even (a,P;E,L) + P CP odd (a,P;E,L) (28) 
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and for antineutrinos 



P{p a ^v p -E,L) = P CPeven (a,P;E,L)-P CPodd (a,P;E,L) (29) 

As an illustration, the probability for to v e conversion in vacuum, assuming three 
family mixing and CP-violation, is given by: 

P C p odd (fi,e) = 2cos(^3) 2 sm5 m sin(2^)sin(^3)sin(2^ 3 ) x (30) 

,AMi 2 2 I, AM?oL AM|oL 

sin — 7^=—) sin ^— ) sin ^— ) 

AE AE AE 

where the CP-phase in matter is: 

e - iSm = (G^-F^)s 23 c 23 + GF(4 3 -sj 3 ) 



'(G>slz + F 2 c| 3 + 2GFc 23 s 23 c s )(G^ 23 + F 2 s| 3 - 2GFc 23 s 23 c s ) 

Here, G and F are given in the Appendix. 

In case of degenerate mass differences, (Am^l ~ |Am| 2 |, the S phase is significantly 
modified from its original value in vacuum for 5 > tt/A radians, while in the case 
lAm^l > lAmlx! the CP phase is almost unaffected by the presence of matter. 

In addition, neutrino propagation in a dense medium makes more difficult to ex- 
perimentally extract a genuine CP violation signal, since the asymmetrical behavior of 
matter, with respect to neutrinos and anti-neutrinos, induces fake CP violation effects. 



10.1 CP statistical significance 

To evaluate the sensitivity to this kind of measurement as a function of the oscillation 
parameters and of the selected baseline, we compare the case where CP violation is 



maximal (5 = 7r/2) and the case of no CP violation (5 = 0). Figure |13| shows the 
sensitivity 

a , r ,,_ N(6 = n/2,E i )-N(6 = 0,E i ) 

b C p{Ei) = == (62) 

/N(5 = 0, E t ) 



i.e. the difference between the two extreme cases, divided by the statistical error (being 
N the number of events in each energy bin, for a given value of 5). We can see that for 
the baseline of 732 km, this quantity is positive for almost the full energy range, so there 
is no real shape variation in the spectrum, and the CP effect is similar to what would 
be obtained with a change in the angle 8i 3 . On the other hand, for larger baselines 
this curve crosses the zero, and the CP violation produces a visible deformation of the 
energy spectrum. 
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Figure 13: Statistical significance of CP violation effects for three possible Am^ values 
mass at the three considered baselines. N(S = 7r/2) is the number of expected events 
assuming maximal CP violation in vacuum and N(S = 0) corresponds to the number 
of expected events in absence of CP violation, for 2 x 10 21 muons 
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10.2 Fitting of the 5 parameter 

To evaluate the precision reachable on the measurement of these three oscillations 
parameters, we perform fits assuming the following reference values: 



Am3 2 


= 3.5, 5, 7 x 1(T 3 eV 2 






Am 2 12 


= 1 x 10" 4 eV 2 






sin 2 8 23 


= 0.5 


#23 


= 45° 


sin 2 2#i3 


= 0.05 


#13 


= 6.5° 


sin 2 6*12 


= 0.5 






5 > 









Figure [14| shows for the baseline L=2900 km, the expected visible energy spectra 
for electron and wrong sing muon class for the two cases 5 = and 5 = tt/2, and 
their difference in terms of number of events. As expected, in absolute values, the CP 
violation affects the two classes in a similar way, but the effect is much more visible 
and significant for the wrong-sign muons, since the total number of events is smaller. 

Figure [14] is normalized to 10 21 /j + decays, since this is the minimum amount of 
decays required to produce a statistically significant observation of CP violation. 

For the largest baseline L = 7400 km, CP violation effects are almost undetectable 
even in case 5 = tt/2 since, wrong sing muon appearance is largely dominated by 
matter effects. At L = 732 km, the effect of a non-vanishing 5 is more striking thanks 
to the higher event rates and the smaller neutrino path in matter. However, as we 
formerly pointed out, this effect is similar to the one produced by a smaller value of 
#13 and therefore, both effects cannot be disentangled with a single measurement at 



this baseline as shown in figure [15], where the correlations between 5 and #13 prevent 
a precise determination of any of them. Nonetheless, a measurement at L = 7400 
km where CP violation effects are negligible and the most accurate determination of 
613 exits, can be combined with data collected at L = 732 km to produce a precise 
determination of S. 

Another possibility to unveil the existence of CP violation is to perform a single 
measurement at L = 2900 km. As shown in figure IT3L at this distance the effect of 



S 7^ is twofold: not only the event rate is modified but also the spectral shape. This 
last effect cannot be produced by a change on 6* 13 . Figure [T^ shows that at this distance 
the correlation between S and 613 has diminished and therefore a better determination 
of the parameters can be achieved with a single measurement. 

For a baseline of L = 2900 km, matter effects are not so strong as for the longest 
baseline. The possibility of measuring the CP-violating phase 5 is not spoiled by the 
fake asymmetries due to the matter interactions, but both effects, even if correlated, 
can be measured at the same time. This is shown in figure [16], where the result of 
a simultaneous fit to the average matter density p and the phase 5 is presented for 
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Figure 14: Energy spectra and differences between no CP violation (5=0) or maximal 
CP violation (5 = tt/2). The two upper plots refer to the electron class, the lower ones 
to the wrong-sign muon class. The plots on the left show the energy spectra for the 
two cases, the ones on the right the different in number of events as a function of the 
energy. The helicity of muons circulating in the ring has been chosen in such a way to 
enhance the effect, i.e. negative muons have been used for the upper plots, and positive 
ones for the lower plots. 
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Figure 15: Correlation between #i 3 and CP phase 5 for two different baselines and 
2 x 10 21 decays. 



29 



S 2.5 
2.2.25 



2 

1.75 
1.5 
1.25 
1 

0.75 
0.5 



21 

L—2900 km, 2 x 10 \1 decays 



Density known at 3% 




J I I I I I I I I ! I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 1_ 



(68% C.L.) 
(90% C.L.) 
(99% C.L.) 

I 



2.25 2.5 2.75 3 3.25 3.5 3.75 4 4.25 4.5 



p (g/cm ) 



=0 2.5 
2 2.25 
10 2 
1.75 
1.5 
1.25 
1 

0.75 




\90im 

C.L.) 











Density free 



0.5 El 

2.25 2.5 2.75 3 3.25 3.5 3.75 4 4.25 4.5 







P (g/cm ) 



Figure 16: Correlation between the average matter density p and the CP phase 5 for 
L=2900 km. In the lower plot we leave Earth's mean density as a free parameter in 
the fit. In the upper plot we assume that density is known within 3%. We see that the 
presence of matter does not spoil the possibility of performing a measurement of S. 
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CP-violation with quasi-elastic events 
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Table 8: Expected QE electron events (iV e / e ) at a baseline of 2900 km in case there 
is no CP violation and in case CP violation phase in vacuum is maximal. The last 
column shows, in number of sigmas, the statistical significance expected for the QE 
electron sample in case CP is violated in the lepton sector. 



the two cases where p is left as a free parameter or it is known beforehand with a 3% 
accuracy. 

In case no CP violation is observed, the allowed 5 values at 90% C.L. as a function 
of Am 2 2 are shown in figure [17] for two different muon normalizations. On the other 
hand, if a significant effect is detected, the precision achievable on the measurement 
of the CP phase 5 is shown in figure [18|. We fit S using all event classes, leaving the 
five parameters governing the oscillation free. Assuming a reference value of 90° for S 
and 10 21 muon decays of each polarity, we get: 90 ± 15°. Therefore a precision around 
20% is expected. Finally, we note that the change on the expected precision for 6>i 3 is 
negligible when Am 2 2 and 612 are included in the fit. 



10.3 Use of quasi-elastic events 

Quasi-elastic events can also be useful to spot the presence of CP violation. Table § 
shows the expected number of QE electron events after kinematics cuts (a back to 
back electron-proton final state with proton kinetic energy in excess of 100 MeV). 
Three different assumptions for the total number of muon decays and the baseline of 
2900 km have been assumed. In case CP is conserved, we expect 35 (96) quasi-elastic 
electron events for Am 2 2 = 3.5(7) x 10 -3 eV 2 and 10 21 muon decays of each polarity. 
In case 5 = tt/2, we expect 26 and 85 events for the two mass differences considered. 
The effect is at the one sigma level. To obtain a statistically conclusive signal for CP 
violation would require more than 10 21 muon decays. 
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Figure 17: 90%C.L. sensitivity on the CP phase 5 as a function of Am^ for two 
different normalizations: solid (dashed) line corresponds to 10 21 (10 20 ) muons decays 
of each polarity. 



32 



,21 



7- 



6- 



5- 



4- 



3- 



2- 



1 - 



E ^=30 GeV, L=2900 km, 2 x 10 \1 decays 



All events 
Only (I 

Electrons and no leptons 



5m 2 23=3.5 x 10" 3 eV 2 5m 2 i2= 1 x 10" 4 eV 2 
sin 2 @ 23 =0.5 sin 2 e i2=0.5 sin 2 2 @ 13 =0.05 5=7t/2 





1.2 1.4 1.6 1.8 

1 Dscan DELTA, all other free 



2 

5 (rad) 



Figure 18: Measurement of the CP violating phase 5 for 2 x 10 21 /i decays and L = 2900 
km. The two mass differences and the three mixing angles have been left free during the 
fit. The reference value is S = 7r/2. The three curves are obtained using respectively 
in the fit all classes, only the muon classes, or only electrons and NC-like events. 
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11 Conclusions 



In this paper, we tried to understand in deeper detail the capabilities an experiment at 
the Neutrino Factory, with the aim of exploiting as much as possible the possibility of 
studying several neutrino transitions at the same time. For this reason, we have con- 
sidered as a baseline detector a large Liquid Argon TPC with external muon identifier, 
the most versatile design proposed so far for large neutrino experiments. 

Assuming an oscillation scenario favoured by present experimental results, the lead- 
ing oscillation would be between the second and third neutrino family, which are maxi- 
mally mixed. Therefore, a very precise determination of the parameters governing this 
transition, # 2 3 and Am| 3 is essential also for the understanding of all other processes. 
This is mainly achieved using the information coming from the disappearance (right- 
sign muon class), provided that the baseline and beam energy are chosen in such a way 
that the first oscillation maximum is visible as a dip in the oscillated spectrum. 

The maximal sensitivity to #13 is achieved for very small background levels, since 
we are looking in this case for small signals; most of the information is coming from 
the clean wrong-sign muon class, and from quasi-elastic events. 

On the other hand, if its value is not too small, for a measurement of #13, the 
signal/background ratio could be not so crucial, and also the other event classes can 
contribute to this measurement. 

Like for a B-Factory, a //-Factory should have among its aims the overconstraining 
of the oscillation pattern, in order to look for unexpected new physics effects. This can 
be achieved in global fits of the parameters, where the unitarity of the mixing matrix 
is not strictly assumed. Using a detector able to identify the r lepton production 
via kinematic means, it is possible to verify the unitarity in — > v T and v e — > v T 
transitions. For this latter, the possibility of a kinematical r identification for wrong- 
sign muon events could allow for the first time a clear identification of this type of 
oscillations. 

The study of CP violation in the lepton system is a very fascinating subject, and 
probably the most ambitious goal of this kind of machines. It will only be possible for 
high beam intensities, and if the parameters governing the solar neutrino deficit are in 
the region usually indicated as large mixing angle MSW solution. Matter effect can 
mimic CP violation; however, a multiparameter fit at the right baseline can allow a 
simultaneous determination of matter and CP-violating parameters. 

Also for CP-violation measurements most of the information would come from the 
wrong-sign muon class, but since in this case the electron class would also be affected, 
the study of these events (and of the very clean quasi-elastic interactions) can provide 
essential cross-checks for these delicate measurements. 
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A Oscillations in matter 



In reference p3| , the authors compute analytic expressions for the mass eigenvalues, 
mixing angles and CP-violation phase in matter, assuming the mixing matrix U is 
parametrized "a la CKM" . Since several missprints were observed, we reproduce here 
the corrected expressions: 

For neutrinos, the mass eigenvalues in matter Mi, M 2 and M 3 are: 

Ml = m \ + -- -VA 2 - 3BS - — VA 2 - WyfTHP (34) 
3 3 3 

Ml = ml + j - ^VA 2 - 3BS + ^-V A 2 - 3BVT zr S 2 (35) 
A 2 



M l = ml + - + -VA 2 - 3BS (36) 



where 



A = Am 2 21 + Am 2 31 + D (37) 

B = Am^Am^D[Am^3 + Am^(44 + 4)] (38) 

C = DAm^Am^44 (39) 

D = 2V2G F N e E for neutrinos (40) 

\l ( 2A 3 - 9AB + 27C\] 

b = cos - arccos = — (41) 

[3 1, 2y /(A* - 35)3 )\ 

and Cij = cos dij and s^- = sin 6^-. The mixing angles and CP phase in matter are: 

• 2 m -(M*-aM 2 + f3)AM 2 31 

12 ~ AMUMf - aM 2 + (3)- AMl(Mi - aM 2 + (3) K ' 

sin 2 ^ - M$^Ml±l (43 ) 

" 13 " AMIAM 2 2 {A6) 

■2nm G 2 s 2 23 + F 2 c\ 3 + 2GFc 23 s 23 c s 

Sill V 23 — Q2 + p2 V^ 1 ' 

e - iSm = (G 2 e-^ - F 2 e^)s 23 c 23 + GF{c 2 3 - s 2 23 ) 

^{G 2 s 2 23 + F 2 c 2 23 + 2GFc 23 s 23 c s )(G 2 c 2 23 + F 2 s 2 23 - 2GFc 23 s 23 c s ) 

where 

2 2, 2/2 2 , 2 \ i 2/ 2 2 , 2 \ { ac\ 

a = m 3 c 1 3 + m 2 (c 1 3C 12 + s 1 3)+m 1 (c 1 3S 12 + s 13 ) (46) 

a 22/ 22, 22\, 2 22 1 ai\ 

P = m 3 c i3l m 2 c i2 + m i s i2> + m 2 m v s 13 (47) 
G = [Am 2 31 (M 2 -m 2 x - Am 2 21 ) - Am 2 21 (M 2 -m[- Am 2 31 )s 2 12 ]c 13 s 13 (48) 

F = (M 2 -ml- Am 2 31 )Am 2 21 c 12 s l2 c 13 (49) 
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For anti-neutrinos, we must replace D by —D. 

In the case of one-mass scale approximation, the oscillation probabilities are simply 
The oscillation probabilities are: 



P(v„ 



v e , E, L) 

E , L ) 
i/ T , E, L) 

E, L ) 
v T ,E,L) 



1 - sin 2 (20 
sin 2 ( 2^3 ; S11! 

cos 



'13J 

sin 2 (2^ 3 ) 
1-4 cos 2 



m ) A 2 

13^32 

23/^32 
2(Qm\ a 2 
23^32 

2/ 



sin 



^3) 



cos 4 (^)sin 2 (2^ 3 )A 2 2 



cos 2 (^™) sin 2 (^ 



(50) 



A 2 

^32 



where Ag 2 



sin 2 ((M| - M$)L/AE). 
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